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IE INTRODUCTION 


Free Electron Lasers (FEL) may have an impact on the 
United States Navy, not just as a weapon, but in many fields 
of science such as Medicine and Physics. The tunability and 
high power output of a FEL makes it useful for a broad range 
of applications. 

As a weapon, the tunability allows the user to select 
the frequency that will propagate most efficiently through 
the atmosphere. Also, the FEL has the potential to produce 
the power level needed for use as a weapon. 

The tunability allows a doctor to select not only the 
wavelength, but also the pulse-repetition rate of the laser 
during surgery. The FEL is able to take advantage of 
resonant modes in a molecule to allow the doctor to 
precisely remove millimeter layers of skin, fat, muscle, 
cornea, or tumor with no collateral damage to surrounding 
tissue. The FEL "scalpel" may revolutionize the operating 
TOOM: 

The FEL is also the most viable path for achieving a 
coherent x-ray source. There are several projects working 
towards a high gain FEL that can achieve 1 angstrom x-rays. 
The uses for coherent x-ray sources are numerous. One use, 
which has a direct effect on the Navy, is to probe the 
reactions that occur in the initial moment of a nuclear 


blast. 


There are numerous other Uses for Ере КЕП qq Almesra every 
field of science is finding new applications for the FEL, 
Ben probing Ligheeineeractioneswitiebiolesıcaleruneenons to 
studying Quantum Electro-dynamics. 

This thesis covers several different FELs which have a 
wide range of applications. The second chapter of this 
thesis covers the background theory that will be used in the 
rest of the thesis[reference (1)]. In Chapter III I 
calculate analytically the height of a tapered separatrix. 
In Chapter IV I analyze various high gain FEL's and use some 
dimensionless parameters to develop a simple model for each. 
The final chapter I simulate the Darmstadt FEL using a 
computer model and analyse the results. This analysis 
consists of simulations at various desynchronisms that 
determine power and gain expected in the FEL. There are 


instabilities that may be experimentally observed as well. 


II. BACKGROUND THEORY 


The physics behind the Free Electron Laser (FEL) is 
based upon both spontaneous and stimulated emission from an 
electron. The two major components of the FEL are the 
relativistic electron beam and a static alternating magnetic 
field. The relativistic beam of electrons is accelerated in 
the transverse direction by the magnetic field. This 
acceleration initially causes spontaneous radiation from the 
electrons. Due to relativistic effects, the emitted light is 
radiated into a cone along the axis of the laser. The key to 
an operational FEL is the interaction between the electrons 
and the light field. The kinetic energy of the electron beam 
is converted to the radiation energy of the optical field. 
The electron pendulum equation and optical field equation 


describes the coupling. 


А. INITIAL CONDITIONS 


A classical approach has proven to be the most fruitful 
and experimental results agree with the classical 
results.[reference 1] 

de Undulator Field 

The magnetic field of the FEL is created by a chain of 
fixed magnets, called an undulator. There are two basic 
types of  undulators: helical and linear. The linear 


undulator is a series of bar magnets with the polarization 


in alternating directions, so that an alternating magnetic 
field is set up. Figure il depiets a linear undulater with 
the z direction ís along the axis of the laser. The magnetic 


field for the linear undulator is given by 


В = [В„, В, в. | = в [o, sin(k,z) cosh(k,y), cos(k,z) sinh(k,y)| , CD WM 
where k, = 2x /A, is the undulator wavenumber. 
mirror 
linear | 5) Å 
undulator a 
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Figure II-1 Linear Undulator (after reference 7) 

The helical undulator consists of two concentric 
currents that create a magnetic field that rotates spatially 
along the length of the undulator. For the helteal undulator 


the magnetic field is given by 


B =>), [cos(k ‚z),sin(k ,z).0] (II-2) 


25 Initial Optical Field 


To begin analyzing FEL operation we need to assume that 
the FEL is already operating and there is already a 
radiation field present. The optical field for the linear 
undulator is assumed to be a plainer, plane parabolic wave 


with the electric and magnetic fields given by 
Ез = Elcos y,0,0] and В. = Е[0, соз w,0] respectively, where 
w = kz - ot +6. The electric field amplitude E is given in 
cgs units with ф being the optical phase and k - 21/ X is 


the optical wave number, 0 - 21C/À is the optical angular 
frequency. 

The initial optical field for the helical undulator is 
assumed to be a circularly polarized plane-wave. Hence the 


electric and magnetic fields are given by 


Es = Elcos W,— Sin y,0] and Bs = E[sin MES wol, respectively. 


В: ELECTRON DYNAMICS 
T The Electron Phase 


Before moving on, we need to define the electron phase. 
The electron phase is the key to coupling the optical field 
with the bunching of the electron beam. In the combined 


undulator and light field the electron phase is defined as 
C(t)=(k+k,)2a(t)+ot. The electron phase definition involves 


two fixed wave numbers k = 2n/i, describing the light wave 


апа k, = 2n / A4, describing the undulator, with the only 


variable being z(t). The electron phase follows the electron 


position. The time spent by the electron in the undulator is 


At=L/ß,c=L/cz10”s, where L is the length of the undulator 


and the speed of the electron (v,=cB,-c) is approximately the 
speed of light. This time is extremely small, so it is 
easier to introduce a dimensionless time т апа а 


dimensionless time derivative 


ВЕ et De 


de = — апа (3) e. ET 





The dimensionless time t ranges from 0 at the beginning of 


ene uUndulatonrmee .) ас сееп л шш е undulator (Osts1) . The 


electron phase velocity is then defined as the derivative of 
the electron phase 


90) È = afk + xp, - x] (11-4) 
Ст 


To study the operation of the FEL, it is most 


ve) = 


informative to plot the phase velocity versus the electron 
phase in what is called phase-space. This will be discussed 
further in subsequent parts of this chapter. 

ze Wavelength of Light in FEL 

The traditional laser requires a bound electron to 
transition to a lower energy level and radiate. In the FEL, 


a relativistic beam of free electrons amplifies the light 


wave. The undulator wavelength A, is Lorentz contracted to 


, 


A, =A,/y in the electron frame. The emitted optical field 


is Doppler shifted is also 


A = (1 -— Ba = 2y. EIS) 
At resonance, there is optimal energy exchange. When 
A, =A in the beam frame. When converted to. the lab. frame 


Бае Теспк ст лоп becomes 


NEUE d 
А = А — (II-6) 
ду“ 


where the dimensionless undulator parameter K is defined as 


еВА eB 
pau — = GT- 


лал kme 





B is the RMS magnetic field, e is the magnitude of the 
electron charge, m is the mass of the electron and c is the 
Speed of light. 


SX Electron Paths 


To analyze the electron trajectories in the FEL, we 
will look at the linear undulator. The magnets are aligned 
along the y direction. An electron that moves off axis in 
the y-direction feels a stronger magnetic field because it 
moves closer to a magnet. The electron is pushed back toward 
the axis creating a natural focusing. The focusing of the 
undulator causes the electrons to wiggle over many undulator 


periods in the focused direction. The wiggle is called a 


betatron oscillation: Toõoanalyze betatron effectes, we will 
consider a perfectly injected electron in a linear undulator 
that focuses in the y direction: Perfect injection means 
that the electron is injected into a perfect sinusoidal 
trajectory with no angular spread. 

The three spatial components are of the Lorentz force 
equation with the magnetic field for linear undulator from 


equation (II-1) as follows 


В, = — еы В. sin(k,z) cosh(k,y) - В, cos(k,z) sinh(k,y)], (II-8) 
ymc 

: eB | 

B, ten В. cos(k,z) sinh(k,y)], i digi 9 

ymc 

: eB | 

В. = - — lp. sin(k,z) cosh(k,y)] . (IT-10) 
ymc 


In the next section, we will derive these equations in 
greater detail for the helical undulator. The p, term has an 


exact solution, and perfect injection implies the constants 


of integration are set equal to zero, so 


ar 


B, = cos(k,z) cosh(k, y) . (II-11) 
y 





Substituting p, into the other two equations leads to 








2 
В, + - СЕ cos'(k,z) sinh(2k,y), (II-12) 
y ү? 
апа 
2 
ß, = На ЕЙ (PE 
y 


Since the y and z components are coupled, we add the 
restriction of perfect injection in the y direction and 


analyze the motion in the xz-plan. Hence, 


J2K 





B, = ш a (ТА) 

ВОКС) АЕ 

Kk 

б. = E — sin(2k,z) (II-16) 
y 


tie neweuxdcst operating FEL Ke andi е hal 





K/y««1. Expand in powers of (K/y), so the x and z components 
become 
— - — 
z(t) - ерес = sin(2k,cB.t] ar or, 
4у 
K AÐ 
x(t) = - —~ sinlk.eB.t) Ps F рео) 
of 


where ß, is the average velocity of the electron beam in the 


undulator. Then converting to dimensionless time T gives 





= Ko e 
ZI Bi = 2 sin(4nNt) +..., CE CEST 
2 Өлү 


Kh, 


x(t) = - os 


From these two equations, we can see that the amplitude of 





sin(2nNt) +... . pM 


the transverse oscillation is KA, / J2ny where the electron 


oscillates once every undulator period. In the z direction, 
the electron oscillates twice every undulator period. 

Having looked at the oscillations in the x and z directions, 
we now go back and assume non-perfect injection to see what 
happens in the y direction. The oscillations in the y 
direction are called betatron oscillations. Assuming small 


sekan сс that ky <Il ага sinh(2k,y) ~ 2k,y we.. Can 


rewrite equation (II-12) and convert to dimensionless 


parameters to get 


se Ae v. 
y(t) 2. - ————— cos" (2nNx)y() (TTE2 I) 


Averaging over several undulator periods, since one betatron 
Eaa 15 50010411) over many undulator periods, mime 


have 


as K’Ik 2 
y(t) = — ——— yí?) = -о; уй). er ze 
T 


Equation (II-22) is the equation of motion for a simple 
harmonie повел Райони The solucion iS a well-studied 


sinusoidal wave, with the dimensionless betatron frequency 
ОЗ БЕ Е 25м оре Mumber of betatron  osciMations along 


the undulator is 


N = —. (ІІ-23) 
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Се рс пе is Nel, Since N=100; K=1, and у-100 аге 
typical values. FELs with larger values will be discussed 


later. 


Ge ELECTRON PENDULUM EQUATION 


The dynamics of the FEL may be described very simply by 
Кс = рс Ыл от. тос ЕЕЕ =Еер л е па 
this equation is to consider an operating FEL where an 
optical field is already present. This means the electron 
experience forces from both the optical field and the 
undulator field. 

The helical field is much simpler to work with 
mathematically and the final formulas can be converted 
easily to the linear case. Hence, we will just use the 
helical magnetic field. 

The forces from the undulator field are governed by the 
Lorentz force equations (cgs units), 

E А sec) = ei Bx), (II-24) 


and the energy equation 


2 
a(yme | - -е ср.Е, (ІІ-25) 
dt 
where 
Е 02 
— = De p^ (II-26) 
y 


me 


The momentum is P=ymv=mßc, its energy is ymc’, and the 
charge magnitude е-е|. We have five equations and four 


unknowns xX, y, 2 апа y. 
Substituting the magnetic field from equation (II-2) 

into cwatıons - (ІТ-24 a25) “апа СБ-26) уе йы the 

following equations 

E. 

abi.) — ЕЙ - D. (cos w,- sin w,0) + 


dt me ; 
p.B(- sin(k,z) cos(k,z)0)] 


(inm 27) 


alyp.) e | 
Z CERE E A 
dt mc | б. со Рат v) 5 TT) 


В, sin(k,z) - B, cos(«.2)) 





d e 
aly) = - — gp, cos v - B, sin v], (11-29) 
sl mc 
where В, = B. + В, and W=kz- ot +6. There are тоге 


equations than unknowns, so we can ignore the p, component. 
The electrons are traveling near the speed of light 


along the z axis, so we can estimate that p,«1. This means 


Einae ans ss 


components of  (II-22) can be ignored. Thus, 
абв.) Ser [p B(- sin(k,z) соз(к,2),0)[, (II-30) 
de mc 


which implies that the static undulator field controls the 


transverse motion of the electron almost entirely. We 
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сс лет Еру snsseesendMEylelding the following 


results 


— 


В, + - 2 [cos(k,z), sin(k,z),0] . 
ii 


SübStituting (ІІ-31) anto (IT-29) yields 


үз cos + 4) 
ymc 


where Č + ọọ = y * kz (k-k,z- ot - $. From (II-31) 


cudecsubstscubig inte to (IT-26) 


So that 


2 


2 


eK) 


~2 


_ 


Á 


N 


Taking the time derivative of (II-35), 


а 


и ва 


From the derivative of the electron phase, 


С = (ко + к„)> — @ 
апа 

er (tee hee 
so that 


(k - k)B- o, 


155 


we have 


(II-31) 


CIT-92) 


CIS 


(II-34) 


(II-35) 


(II-36) 


(II-37) 


(II-38) 


hee 
5 = 
ee тке (II-39) 


(к + кје 


ОШО Еп тпс тпс От -зо we gert 


ки ee IP: | (II-40) 
ү 1+к (1 + кк + к.) 


Now equating (II-32) апа (11-40), we get 





2 .. EK 
ше 


(i+ Kk # ke 
and solving for ee we have 


- eos uuu (II-42) 
Un D. yme 


sance k>>k, we can use the approximation (k+k,)=k. From 
Sct Rene = CK, = ск + K?) / 2y?, and using p,«1, 


equation (II-42) reduces to 


s — 20,eEK 
C = S adds cos (6 + ф) писа 
YE 


Converting to dimensionless time with dt-Ldt/c, and defining 


the dimensionless optical field magnitude as 


AnNeKLE 

pc (II-44) 
y me“ 

The pendulum equation becomes 

б = lalcos(& * 9. (II-45) 


The pendulum equation (II-45) is the same form for a linear 


undulator. However, for the linear undulator, the optical 
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field la| is modified by substituting K[J,(E)-J,(E)] for K, 


where E-K^/2(1«K^) and J,(&) and J,(&) are Bessel functions. 


De OPTICAL WAVE EQUATION 


The pendulum equation describes how individual 
electrons are affected by the laser field. To describe how 
the laser field is affected bv the beam of.,electrons, we 
next derive the optical wave equation. Begin with the 


inhomogeneous wave equation for the optical vector potential 





F 

мо до (11-46) 
с ас c 7 

where 

B 1 OA | 

о ee E[cos ШЕП у,0], (II-47) 
CUTE 

pe VAS [sin U, COS wo], (II-48) 


y = kz- ot +0. When the amplitude and phase of the laser 


beam change slowly with time, we can write the vector 


potential in the form 


— 


Ze £ Ez, t) [cos (y) ,- sin (y) ,0] . (II-49) 
W 


The form in equation (II-49) assumes that the electron beam 


is small compared to the waist of the optical mode, so that 


the transverse direction is not important. A is dependent 


only on z and t. The spatial derivative is 


5 


CD 
№ 
Зе 





Q 
O 
(Л 


(y) ,- sin (wy) ,O] + 


г 


С) 
N 
| 


= 


CD 

Er] 

To 
+ 

2 





2 
í [sin( (WIT COS) „0) + 
да 
ОЕ Ó 
— [ + 2 [cos (Ч) ,- sin (vw) 0] + | Crt sie) 
02 02 
БЕ 








[соз W Sinay) „| - 


8 2 
= [- ӘНШІ Cos (1) 0) 


mE cb rie xp ge 
EU d 
n г 


and the time derivative is 


), cos (w) ,0] 





— o [cos WV Siny OS 


ad“ 
ot 





--- - o [ces (W) Sn) „0) : [cos CU enm T 0] E - 515 


- E — o| Sam) eos) 0) 


From the assumption that phase and amplitude varies 


slowly, ðE/ðz << КЕ, 02/02 << ko, OF/Ot << WE, and 06/dt 


<<00. From these approximations, we find that 





= mo ОЕ | ole 
i X == a Еа eost), sin(y), 0] (II-52) 
ZI AL sin(w),—cos(y), 0]--227,. 
07 с 


To further simplify, we use a coordinate transformation, 


where the position z moves along with the optical field at 
the speed of light c, so that 2 = Z + ct. Also, we Convert 


to dimensionless time, t=ct/L to get 
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Ó О О 2, 


62 côt) côt La 
and 
І ОЕ lo 4л- 
[сому] ds sin(w),-cos(w),0|= тола ; (II-54) 
L or L ox C 


Each electron adds to the total current, so the current 


density contribution from the i^ electron is summed to get 


1, =-ecS 8,8 x - x) (D SS) 


where the r is the trajectory of the i" electron and 8° 
Ше а three dimensional” Dirac Delta Finetzonee ubseirutiıng, 
the transverse velocity By = -K / y(cos(k,z), sin(k,z),0) into the 


current gives 





ү. Jie» - ш ре (kaz); sin (k,z)] 59%; = г). (11-56) 
С і i 


Therefore, the wave equation is 


HZ (ume eser UA „ol + 22 Sumus eS „ol - 
От дт 


(II-57) 





Дан one [cos (k;z) , sin (k,z)] 59%; u =) 
1 ү 


Both sides of (II-57) averaged over a volume element а\,, 
which is smaller than the coherence volume, yet larger than 
a wavelength of the light. The current density over this 


small volume element is constant, so the summation over 1i 


electrons >» (- 2neKL / у)соз (оаа и (k,z)] 60) can be replaced 


ша 


with a charge density p times the average of sample 


electrons in the beam, the average is symbolized by F So 


now we have two uncoupled equations 


СЕ 2nepLK 
л C iU) (II-58) 
Ct | 


Ó 2 LK 

с И 0) (II-59) 
От 
where U-k,z-$-9. Defining the dimensionless current 


j = 8n’e’pNK’L’/y’mc’, and recalling the dimensionless optical 


тема lal = 4nNeKLE/y’mc, the change in electric field 
amplitude becomes 


9E . 8 








- -Хсоз (С + ф)). (mio) 
OT 
Since the electric field is complex, then a = ае“, and we 





arrive at the optical wave equation 


a--je^). (ISP 
Equation (II-61) is derived for a helical undulator. 
However the same form of the equation applies to a linear 


undulator provided K is replaced by K{[J,(&)-J,(&)]. 


Е. GAIN MECHANISMS 


For the FEL to be operational, the there must be an 
increase in power of the optical field each pass through the 


undulator. There are two regions of optical gain: the low 


ща 


jaime меса елт) -and tbe high gain region ()>>я). To 
understand the mechanisms behind gain we first develop the 
concept of an electron phase space plot. 


та Phase Space 


The electron phase space is a plot of an electron's 
phase velocity, v(t) = á = (к + k,B, - k] and the electron 


phase €(t)=(k+k,)z(t)+@t. In the low gain case а * 0, which 


means that equation (II-61) evolves slowly so that lal and 9 


are approximately constant. Multiplying both sides of the 


electron pendulum equation (II-45) by € yields 


a = а| соз б + ф)ё = vv (617262) 
(OUS 

vv-lalcos(& + $94 = 0 ; MENS) 
so that 

$T - lasing + o) - 0. (II-64) 
ar? 


This means that v’ - 2lalsin(& + $) = constant, and 
vo =v, ~ Jalkin + 6) - sin(C,)}. (II-65) 
Where v,=v(0) and €,=C(0) are the initial phase velocity and 


phases. Equation (II-65) is а statement of energy 


conservation in the dimensionless coordinates € and v. The 


T3 


equation also identifies paths in the phase space, víð), 
that describes the electron motion. In Figure II-2, several 
phase space paths are plotted. 6 refers to taper and will be 


discussed in the next section. 


FEL PHASE SPACE PATHS... 


а=1 ф=0 б=0 


dha Open Orbits a" 





4Уа 





-д/2 C 31/2 


Figure II-2 Phase Space Diagram 


2. Separatrix 
Electrons stay on v(€) paths in the C€C-v plane. These 
paths are based upon the initial phase velocity wv, and 


initial phase Ġa Oof the electron and спе pendulum equation. 
This same equation describes an undamped, undriven pendulum 


where € is the tilt angle and v is the angular velocity. The 
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батса отне есе paths "occur at v=0 “and бел/2#лп 


where n=1,2,3. The stable fixed points occur when б=л/2+2лп. 
In Figure II-2 the stable fixed point is labeled with an 
"S", The unstable fixed points are at v=0 and 6=3n/2#2nn. In 
Figure II-2, the unstable fixed points are labeled with a 
"Jj". At these fixed points, the electron remains stationary 
in phase space. At (v,6)=(n/2,0), the pendulum is hanging at 


| 2/2 


the bottom of the arc. When Ivi»2la , the pendulum has 


enough angular momentum to swing continually over the top 
around in a circle. In phase space, we can draw a separatrix 
which separates the lower closed orbits from the over-the- 


top open orbits. The equation for the separatrix is 

у, 2 2la[1 -sin(C, « 9)]. тт 
In Figure II-2, the separatrix separates the open orbits, 
shown as a gray area, from the closed orbits, shown as a 


white area. 


3% Bunching 


Consider a collection of electrons that are initially 


uniformly distributed in phase evolving along the undulator. 


The pendulum equation shows that for -1/2 < G < л/2, G»50 and 


the electron phase accelerates, while for 7/2 < б < 3n/2, 


C<0 and the phase decelerates. The result is that the 
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electrons tend to bunch about л/2 as they evolve along the 


асо еза 


4. Changes in Phase Velocity 

From the definition Of phase velocity 
v(t) = Lk + k,B, k], we find the change in beam energy due 
to a change in phase velocity. Taking the derivative of both 


Sides near resonance where v=0, k>>k, and N=L/A,, we get 


Ov = LkoB, = NA,Köß,. From equation (II-6) and the resonance 
eondarienm AF: AQ K yay? , we find that 


бу = Ann = on D. (II-67) 


Y 
Thus, an increase in v results in an increase of the 
electron energy, ymc/, and the optical wavelength A. 
5. Gain Formula 
Gain is defined as AP/P, where AP is the change in optical 


power in one pass through the undulator, and P is the power 
at the beginning of the undulator. From (II-67) the average 


change in energy for an electron in the undulator is 

Ayme? є ymc?((v) — vg) / 47N, (II-68) 
where (v) is the average final energy. The number of 
electrons in a volume dV of the optical wave is given by 
pFdV where p is the electron density and F is the filling 


factor. The filling factor F is the cross-sectional area of 
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the electron beam divided by the cross-sectional area of the 


ght beam. The radiacion energy пп а volume dV is given by 


2E’AV/8n in cgs units. Invoking the definition of the 
dimensionless field a and the dimensionless current j, we 
arrive at 


- leravime (o - vam] zi _ у 


G 2 2 
2Е-ау / 8л Ay 


(II-69) 


6. Low Gain 


For low gain, the optical field changes vary slowly so 


we can estimate that a «0 with the initial conditions 
v(0)=v, al0)=a, and G(0)sG,. From perturbation theory, we 


[uude that 


UN VEL 2 [sin(C, * ovs) - sin(é,)| + 
0 
T 
ao u - со5(26, ar Фут) т COSME at cos(2v,t) - (І 0) 


3 
Фо 11 - утз1 (С) соз(6 % Уот) 


Taking take the time average of the phase velocity, we find 
a 2 

oS + та [2 cos(v,t) = и sin(v,t)| T sso Те (como) 
0 


The gain is therefore 





С = 2E (v = v) = ы [2 = 2 соз(у,т) = Vet sin(v,t)]. 2) 


= 


а, 0 


© 
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е. High Gain 
Minne high салш regrmswWE-he Current is very large, 
Jeon. Unlike tne low gain segame, we cannoewignore the rapid 


evolution of the optical field over a single pass. Combining 


the expanded pendulum and optical equation ап тпесеставиеит 


over all the initial electron phases, (ac, , we get 


NO = = Кз: - т!) (ТТ- 73) 


where F(v)- ааа“ ШОУ ‘characteristic function of the 
du cmn is the distribution of the initial 
electron phases v;2-v,«q about v, and [daf(a) 21 : 

The low current gain formula can be obtained from this 
equation by considering a perfect beam where f(q)só(q) and 


ат-т) “а, 11 which case a, can be moved outside the 


integrand. Further integration results in the low gain 
formula shown earlier. 


To find the solution to equation (II-73), we need to 


consider a perfect beam with F(t')= 1 and a resonant optical 


field v,=0. Then, taking two derivatives, we get 
ооо 17 
alt) = > a(t). (II-74) 


lucu scehowuWwPHact ehe Solution CO (II=74) TS 
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a : RT + > 0053143099 ЕО 172043. 
a(1) = à poro (i V3) /2 Е е ЭУ 43) / 2 ET e0/ 9 | (II-75) 


Since for high gain jJ>>n the dominate exponential terms аге 


the fastest growing terms. The following approximations can 


be made for the optical field amplitude, phase and gain: 


alt) = 2 Bu oue. (II-76) 


5 1/3 
Е p 
«9» [3| >" ттш 


сб) оо (TIS) 


The gain grows exponentially along T with a growth rate 


Корок опа to q | 


575 BEAM QUALITY 


Previously, when considering the interaction of the 
electron and optical pulse we assumed perfect injection. In 
reality there are errors involved. The two factors that have 
the most influence are electrons injected off axis at the 


initial position, y,, and electrons injected at an initial 


апе ло Prom Reference "(2)PhecdecreaseSun phosegsue las 


TS 


= 2nN пати De 
Av = —— |K*k * vð. (11-79) 
1 + к? ( 0 Yo Y y ) 
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The change in phase velocity from these errors is made more 
severe with more undulator periods N. A longer FEL is more 
sensitive to injection errors. 


A measure of the quality of the electron beam is the 


smita. rs the кешт шал рост оп 


spread of electrons along у, and 8, is the rms initial 


angular spread of the electrons in the y direction. The beam 


oe oom с EO minimize either yer 8) олет 


sma fixed. The minimum change ang phase 


velocity occurs when the beam is "matched" so that both 


terms in (11-79) contribute equally. Matching occurs when 


Kk,y, « y9,. For a matched beam Av * 2nN,Kke, / y. 


To show the effects of beam quality we need to look at 
the FEL integral equation (11-73). The function f(q) in 
equation (II-73) has a distribution based only on the 
initial distribution of the electrons as they enter the 
undulator. Two examples of distributions are those from the 
electron phase velocity and injection angle. The spread in 


phase velocity can be considered a Gaussian distribution 


about v, with a standard deviation of 6; = 4nNAy/y. The 
angular spread leads to a exponential distribution with a 
standard deviation given by 0, - 4nNy' 0, /(1 + kê). The 


distributions are 


26 





e (II-80 
i 4210, : 
and 
ga’ 
(а) = : (LI=-3% 
Og 
The characteristic function for these two cases are 
oa (TI 32) 
and 
1 
DS —_— (1783) 
Lc dT 


For perfect injection Og=0 and 09=0, we have F(t)=1. When the 
beam is not perfect the characteristic function | F(t) | 


decays with a characteristic 4іпе 1/6; ос 1/69. F(t) depends 


on the shape of the function  f(q), and different 


distributions of f(q) effect the FEL's ability to bunch. 


G. SHORT PULSE EFFECTS 

In the low gain FEL, the modes of light in the 
oscillator can develop independently. Since the electrons 
are relativistic, they travel to, but less that, the speed 
of light. At resonance, one wavelength of light passes over 
an electron as the electron travels through one undulator 
period. Over the total length of the undulator the electron 


slides behind the light by NÀ, called the slip distance. In 


the transverse dimension, the optical mode is roughly the 
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radius of the optical beam which is (NU m These modes 


evolve as determined by the pendulum equation and wave 
equation. 
ІШІ ексек toto evolution ОЕ ИСЛЕ о ЕЕС Е: Шав 


we must set up multiple sites in z along the complex wave 
envelope. The transformation to spatial modes, a-a(z), is 
the equivalent to the transformation to wave numbers, 
a—a(k). The longitudinal position on the z-axis 2 is 
transformed to dimensionless parameters by dividing by the 
slippage distance NA, therefore zZ/N\.  Througout the 


remainder of this paper z is dimensionless parameter. 


The length of the electron pulse 1, becomes the 


e 


dimensionless parameter o0,=1,/N\. The pendulum and wave 


equations of motion with slippage included become 








к= о i) (тте, 
апа 
A ое. (II-85) 


The electron current is j,=8n’e’p(z)NK’L’/y’mc’, where the 


electron density, p(z), is the density at the site z. A 


site in the electron beam interacts with and amplifies N 
optical wavelengths, while the optical wave passes over and 
attempts to bunch electrons in a slip distance along the 


electron beam. 
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chort oul Ses of spontaneous radiation drive a typical 
ГЕ ес аса леп ые length of ithe electron pulse, E 
is comparable to the slip distance, short pulse effects 
influence the interaction. As the electron and optical 
ӘН нес раса over eachother the pendulum equation (11-85) 


sees a continually changing optical field, a, while the 


wave equation (II-86) sees a continually changing current 


аетпе БЕЗИ Чите: пао distribution CE. n In the 


resonator, the optical pulse bounces between two mirrors. 
The distance between the mirrors is S which is necessarily 
greater than the length of the undulator, L. The time 
interval for the light-pulse make one round trip is 2S/c. 
Desynchronism, d, is the difference between the arrival 


times of the electron and optical pulses at the beginning of 


the pass, again normalized by the slip distance NA. If d=0 
at t-0, the electron pulse arrives simultaneously with the 
light-pulse. Moving one of the mirrors slightly by AS 


changes d so that d=-AS/NA. 


Hc TAPER 
I. Purpose of Taper 


In a successful FEL design, the optical field will 
eventually saturate due to over-bunching of the electrons 


along the undulator, thereby changing the dynamics of the 
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electron interaction. As the electron loses energy to the 
optical зела Еве внасе ие осетр шпоуесжәГЕ Шесспапсе Шыр 
Avz4nNAy/yz-2n. One way to extend saturation to higher 
powers is to taper the undulator field. Instead of a 
periodic sinusoidal magnetic field with constant amplitude, 
the amplitude of the field is either increased or decreased 
progressively along the -length of. the -undulator. Either 
shortening the undulator period along the undulator, or 
changing the strength of the magnetic field along the 
undulator tapers the undulator. In either case, acceleration 


of the phase of the electron results. 


This acceleration is denoted by 6, or the taper 


parameter, so that v=6, and 


KAK 
6 = —4xN ——— , (II-86) 
1+ K* 


where the undulator field is tapered by AK over the 


undulator length. The electron pendulum equation becomes 


C=v=5+lalcos( +9). (II-87) 
While the optical wave equation (II-86) remains unchanged. 
25 Taper Phase Space 
In the limit of low gain the phase space paths for the 


tapered undulator are given by 


1 


у = lv," - 25(6 - 6,)- да вт (С + 6) — за (6. + УЕ (II-88) 
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тышеч ые лп for both positive (60»-0)Land negative 
taper (0<0) in phase space is at 

(Со) - (со (-З/а)- 6,0) (11-89) 
The unstable fixed points depends upon the direction of the 
taper. The unstable fixed point for 6>0 in phase-space is 
located at 

(Ev) = (2n-cos™ (= /|a|)- 6.0) (II-90) 
and for 6<0 it is at 

(C,v)= = (- cos '(-8/lal)—¢ 0). (II-91) 
The formula for the separatrix is 

v, = 28(6. - Go) lallsin(G, - 9) sin(G, +9), sn 
where &,=2n-cos"(-ö/a) for ö>0 and &G,--cos'(-8/a) for 6<0. 


Figure II-3 shows phase space plots for negative and postive 
tapers. "U" refers to the unstable fixed point, and "S" 


refers to the stable fixed point. 


Se 


FEL PHASE SPACE PATHS | FEL PHASE SPACE PATHS 


| az20 o=0 бш-2.5Л | а-20 0-0 22.54 
ТЫ ТЕО р = 





-1/2 C 31/2 -1/2 C 37/2 


Figure II-3 Phase Space Plots for Negative and Positive 
Taper 
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III. TAPER SEPARATRIX HEIGHT 
To find the height of the separatrix for the tapered 
undulator it is convenient to make a few variable changes 


and assume that there is no initial phase (which only shifts 


the Separatrix along x) к=, у = у, ф = 0, S 
Co = Xy. Hence, equation (II-95) becomes. 
y? = 26(x - x) * 2la(sin x - sin xj) Сета) 





-1/2 X 31/2 


Figure III-1 Separatrix for negative taper with qz0.4 


Figure III-1 shows the separatrix for a taper 0--2.5 and an 


optical field a,=20. The height of interest is the maximum 
height labeled h in Figure III-1. To find the maximum value, 
we need to take the derivative of equation (III-1) which 
gives 


2yy' = 28 + Jal cosx, suu 2) 
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where y = dy/dx, and set the slope y’ equal to zero 


5 + lal cos x 
V = Ec 3 = 0. (III-3) 


The maximum value for x can be denoted as x and is given by 
Ò + а eos = 0 or sg) (III-4) 
where g = 0 Да! . Note that q is restricted to values between 
Lane). 
A. POSITIVE TAPER 

КО ОШО x. 2r-eesr a) and substituring x for x Gives 


у^ = 28\cos™ (— q) - (2x - сов (- 3)) 4 


Jalsin(cos( a) - sinr - cos" а) 111-5) 


Using the identity sin(?r-x)=-sin(x) we can simplify equation 
(III-5) to get 

y? = 25| cos? - a) - 2л| + Aale sin(cos'(- а))|. (III-6) 
Factoring out |а| we have 

у? = 4alalcos"(- a) - л) + sin(cos?(- a)) . (III-7) 
Using the identities cos (- x) - n 2» - cos (х) and 


sin(cos"(- a)) - J1-q' , we can consolidate equation (III-7) 


LO 


у? = jal- q соз (q) + 4/1 – а? | | (III-8) 
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Taking the square root of both sides, we get the maximum 


value of y to be 


NE ofall a cos" (a) * 41 — ET (III-9) 


ОСИ the separati i h 55 twice the value otf v 


therefore, 
1 
h = aill- q соз (а) + 41 - а | ши i БЕРСІНШІ 


B: NEGATIVE TAPER 


For negative q, x, =-cos (-q). Therefore equation 


(III-1) becomes 


y^ » 26|cos "(7 q) - E cos (= а)) + 


salsinlos"( a) - sin( сов ај. Бан 


Utilizing the identity sin(- x) = -sin(x), it becomes 

у^ = 48|cos" C q)] + 2lal- sin(cos"(- а)) - sin(cos^(- а)) (ІІІ-12) 
Consolidating and factoring out а! gives 

y = dlalla(cos*(- а) E sin(cos"(- а) , (ТЕЕ Э) 


Take the square root of both sides gives 


: 
у = Бы cos (- q) * 41 - qi (Еа) 


and then the height is 


E 
h = á fella cos (- а) + 1 — en (III-15) 
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Equations (III-10) and (III-15) are equivalent with the lal. 


Therefore we can say for all values of q 


п = fall [el cosa) + Mi = | (111-16 


ho | (P 


C. POSTIVE AND NEGATIVE TAPER 
Since the separtrix height for the non-tapered case 


(620) is 


2 
h - a ТТ) 


Figure III-2 shows a graph of equation (III-18). 





Figure III-2 Change in Separatrix Height ó/a 
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I» LINEAR APPROXIMATION 


A Taylor expansion inside the brackets of equation 


(III-18) and using the expansion of arccosine 
cos а) + ——-aq +... (TILI- IB) 


gives 


h = sal - (5 = ] а - Y ns 


which simplifies to 


VEN 


> 
И 


«fe: = 54 + TH (III-20) 


«fa - Z a) = «ai: - 0.785 5 Е ав) 


а 


D 
Il 


From equation (III-21) it can be shown that, as ô — 0 the 


separatrix height һ — 4„а, or the size becomes the same as 


the untapered case. As the optical field a > oo and the 
height h = the same as д-0. Ав «һе optical field 


becomes zero a => 0 then the height h > 0. 
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IV. PARAMETERS FOR USE IN SIMPLE MODEL OF A HIGH GAIN FEL 


To develop a simple model for a high gain FEL we need 
to first take a close look at the parameters that dominate 
the electron interaction with the optical field in the high 
gain regime. Table IV-1 is a table of several proposed and 


well-documented high gain FELs. 


[isses | was [resa mures | mee [meer | mer — 
кє [мез [зо (3500 | 33 | 35 | $5 | 23 
CT ames | 3400 | sooo | soo | eso | 1005 | эзо_ 
am þa | o | 3s | 3 | « [as 
ST C27205 | 35000 | 3305 | 5060 | 3660 | 3000 
x [ = | oas [ocos | oos | 0.08 | 0.33 | 036 
ре ре ре Ги Тежат 
Dom [3 [5 [9*9] 93 
ca [- oe [ame | 3» | 3e | 4 | 
Cw] - | 9 os þa ou [$25 
al Tee sr fe fe | 32 | зе 
[= [тов [ш [зо [зз | 300 | 90 
m| e so ЕЛГЕН БЕГЕН ШЕТЕН БЕЛІН БЕСІН 
eris = 306 | 3006 | ae | ee ШЕТЕШЕЛГЕ 
Се en ооо оо БЕГЕН БЕГЕН ШЕННЕН 
D p em m [ewe[ 1s | x3 | xs | 95. 
ret - [2 jo [eee] es | 0.28 | 


Table IV-1 Table of High Gain FEL Parameters 
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The following chapter is a description and analysis of 
each FEL ain order to characterize them with our simple 


model. 


A. DESCRIPTION OF LCLS AT SLAC 


The Linac Coherent Light Source (LCLS) proposal is to 
use the last kilometer of Stanford Linear Accelerator Center 
Ku ENS. ERES I.E ros A 52-520 
wavelength. An electron gun will inject an electron bunch 
that will be accelerated along the last kilometer of the 


SLAC accelerator to 14.35 GeV energy with a peak current of 
3400 A. The normalized emittance is £,of 1.5 mm-rad. After 


the  undulator, the electron beam will be dumped. The 


accelerator is tunable from ваша EMEN GeV which 


corresponds to x-ray wavelengths of 1.5-15 A. The electron 


=з 


þsam particle densiky 15 р-2 5Х107 сп. іп a beam radius of 


r.=30u. The undulator will be a linearly polarized with a 


Fnit Msg» "meeeriod! A.—3 09cm Memumber Ser periods 


N=3328, peak magnetic field strength B=1.32T for a resultant 
undulator parameter of K=2.6, and a dimensionless current of 
j=27,000. 

The FEL will be a single-pass amplifier starting from 
noise or spontaneous emission; a process called Stimulated 
Amplification from Spontaneous Emission (SASE). Due to the 


lack of mirrors at this wavelength ,numerous proposals are 
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being investigated to guide the light once it leaves the 


PERS 


B. TESLA FEL AT DESY 


The Deutsches Elektronen-Synchrotron (DESY) in Hamburg, 
Germany, is also planning to develop an x-ray FEL using the 
TeV Energy Superconducting Linear Accelerator (TESLA). The 
electron beam will Бе supplied by two 16 km long 
superconducting linear accelerators which were initially 
designed to collide electrons and positrons. 

The specifications for the TELSA FEL are still under 
consideration and the project is being planned 31321 
conjunction with the LCLS with a later timeline. Many 
parameters will be based on experiments conducted in 
preparation for LCLS. The projections are based on beam 


energy of 25GeV and a peak current of 5kA. 


С. ANALYSIS OF SLAC AND TESLA PARAMETERS 
1. Wavelength 


The purpose of these two lasers is to create a source 
of high power, coherent, hard x-rays at a 1 À wavelength for 
numerous experimental uses. The lack of normal incidence 
optics at this wavelength creates many unique problems. Only 
grazing-incidence optics are possible, and they are still 


under development. The limited optics precludes the use of 
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an oscillator contfigurationwand requires wheXFEL to operate 
as a Single pass amplifier. 

2: Reducing j 

The dimensionless current j for both x-ray lasers is a 
large number. There are several factors, which lower the 


effective j for actual use in our simulations. 


a) Fillıng Factor 

In a typical FEL, the radius of the electron beam 
is smaller than that the optical mode. The filling factor is 
the cross sectional area of the electron beam divided by the 
cross sectional area of the optical beam. For LCLS and TESLA 
the cross sectional area of the electron beam is comparable 
soli ross section area of the optical beam, so F= ror 


Kan пе гал а сс ее eron beam r. IS about Oss mm 


30 um so the normalized pulse width 0,*«0.45. The radius of 


the optical field is approximately w,z0.02 m- 2 um. 


Gaussian 
Opical Beam 





Figure IV-1 Gaussian X-ray Beam and Top-Hat Electron Beam 
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b) Betatron Oscillation 


The number of electron betatron oscillations along 
the undulator in the x-ray FELs is small. The number of 
betatron oscillations along the undulator is inversely 
proportional to the energy of the electron beam. If there 
were a large number of betatron oscillations, then the 
average amount of time electron beam spend within. the. 
optical field is of little concern. Betatron oscillations 
cause the electrons to pass from the outside of the optical 
mode, where the fields are weak to the center of the mode, 
where the fields are strong. In order to account for the 
reduced bunching drive of the electrons due to betatron 
oscillations, due to the dimensionless current  j ls 
decreased by half. 


За Rayleigh Length 


The waist size of the x-ray beam is very small with a 
w,220 um radius. A small Rayleigh waist means a small 


Rayleigh length and large diffraction pattern a large 
distance away. The x-ray beam diffracts away from the 


electron beam contributing to reduced coupling. The Rayleigh 
length, 2z = Ты Ади for the LCLS this is approximately 


30m, about one-third the length of the undulator. In a high- 
gain FEL such as LCLS the gain is exponential along the 
undulator. For LCLS, the gain in a Rayleigh length is about 


100m as shown in Figure IV-2. 
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Optical beam spreading 
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Figure IV-2 Discrete comparison of gain versus diffraction 


TO illustrate the reduced role of diffraction of the x- 
ray beam, divide the undulator into three discrete pieces. 
At the beginning of the undulator the optical power is taken 
to be 1 as shown in Figure IV-2 on the left. In one Raleigh 
БӘНИ” тепе eget, Or the undulator, ene саш тє 
100, but only within the region of the electron beam. The 
power in the center is 100, while the x-ray beam outside the 
electron beam has had no gain and remains power 1. The 
strength of the diffracted x-ray beam is only 1% of the 
power of the x-ray beam in the center. In the second 
Rayleigh length, the x-ray within the electron beam region 
has been amplified another factor of 100. The power in the 
center x-ray beam after two Raleigh lengths is 100^. The 


power of the diffracted beam was amplified to 100 in the 
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first Rayleigh length and remains 100. Again this only 1$ of 
amplified x-ray beam within the electron beam region. In the 
third length, the power at the end of the undulator has 
become 100°, while the diffracted light has only power of 
100°. The power of the diffracted region is again only 1% of 
the power of the beam in the center at the end of the 
undulator. 

In reality this is a continuous process along the 
undulator. But, just as for this discrete example, the 
continuous gain within the central area of the electron beam 
far out-weighs the gain at the edges of the electron beam. 
TO account for this reduction in coupling, the value of j 


can be reduced by half, so that considering betatron motion 

and diffraction the final value is  j«6700. This e 

comparable to the other long wavelength FELs in Table IV-1. 
4. Emittance 


The emittance for the high gain FEL can be extremely 


important. The x-ray field growth is exponential with an e- 
folding time оғ т, » 2/G/2)""V3 = 0915 For j-5000. For LSes 
the emittance of 0g=0.7, leads to a e-folding decay time of 


Tt, ~ 1/o, ~ 1.4. Since the decay time caused by emittance 


б 


is greater than the growth e-folding time driven by the 


electron beam, т. = 1.4 > т, Ж 0.1. The net growth rate of 


с 


the x-ray beam is significantly reduced along the undulator. 
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5% Slip 
The electron pulse length is about 186 time longer than 
the slip distance NÀ. The extremely short x-ray wavelength 


makes the slippage distance very small and short pulse 
Ss Tare ENON CXISCent ЕСП РЕ О ое а pulse 
experiences gain, which is proportional to the local density 
of electrons. 

6. Initial Optical Field 

Since ше other laser can radiate near x-ray 
wavelengths, the x-ray FEL must seed itself. The initial x- 
ray field a, for an x-ray laser is necessarily zero, and the 
laser depends on the SASE (stimulated amplification of 
spontaneous emission) to begin. A simple simulation model is 
to start with a very small initial optical field that is 
representative of the spontaneous emission created in the 
Desgsnnind Secemem of thegindulator. Іп the first Soemeron of 
the undulator, the model assumes that an initial optical 
field starts from random fluctuations in the electron beam 
density. The x-ray mode in SASE is known to be linear, 
Starting from random noise, and growing through mode 
сешресісіоі СО спе desired wavelength. Since the field 
growth in the beginning of the undulator is linear, instead 
sigrar along the undulacor, “а 


slightly longer undulator is required. 
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Dr ELECTRON LASER FACILITY (ELF) 


ThegEBectroneLbaser Facility ELE) is an offshoot of the 
Experimental Test Accelerator (ETA) located at Lawrence 
Livermore National Laboratory. The ELF was built to utilize 
the electron beam accelerated by ETA. In each year, 1985, 
1986 and 1987 a paper was published updating that years 
experiments at ELF (references 4,5, and 6). The main 
components of the laser remain essentially the same with 
several modifications from year to year. Here, the three 
phases of the ELF are referred to as ELF85, ELF86, ELF87. 
ELF85 had beam energy of 3.3 MeV and ELF86 and ELF87 had a 
beam energy of 3.5 MeV. ELF85 and ELF86 used a 3m long 
undulator which was extended in ELF87 to 4 meters. The 
undulator for all three is composed of specially shaped 


solenoids with a peak magnetic field of 5 kG. The wiggler 
period remained Á=9.8 cm for all three. The number of 
periods was N=30 and Nz40 respectively. The peak current 
initially was 500 A for ELF85, then increased to 850 A for 
ELF86, and finally to 1000 A for ELF87. The wavelengths 
produced were A=0.76 cm, 0.72 cm, and 0.175 cm. For each 
experiment, the optical field was contained in a 2.9 cm by 
9.8 cm rectangular waveguide and the FELS were operated as 


Single pass amplifiers. 
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E. ELF85 
ШЕСЕ has a dimensionless current ELE The 
тпепотоптезв”реантет ге 0, = АБИ che vns NND eal 


field of a,212 was produce by a 60kW pulsed magnetron. 


Hu Electron Beam to Optical Field Ratio the Filling 
Factor 


A sketch of the cross section of electron beam and 
optical mode is shown in Figure IV-3 and the filling factor 


is about F=0.03. 





Figure IV-3 ELF 85 optical field(white), electron beam(gray) 
diagram 
The electron beam is small enough so that electrons see only 
the peak of the optical beam. There are a large number of 
betatron oscillations along the undulator, unlike LCLS, but 
the electrons continue to see the peak optical field. The 
betatron oscillations increase little effect, since they 
remain within the optical field. There is some radiation in 
the harmonics due to the oscillation, but not enough to be 


significant. 
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2n Slip 


шпшенешесетоп ришестишен20 с пррасе distances in length 
tor EbRssr This icme maller Ehan for the LCLS, but still 
much larger than 1. Since the slip is larger than 1, short- 
pulse effects can still be neglected. 


3% Emittance 
The emittance for the beam was quite large £20.37 cm- 
rad. This large emittance, assuming a matched beam, means a 


090-22, and therefore a characteristic decay time of 


Ts=1/09=0.045. Gain degradation due to the emittance is quite 


Sum fue 


E: ELF 86 

The major difference between the ELF85 and ELF86 is the 
modification in the Experimental Test Accelerator (ETA) to 
use a field emission cathode. This increases the peak 
current and makes the pulse longer and the emittance 


smaller. The undulator dimensions remained the same. 
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T4 Electron Beam to Optical Field Ratio, the Filling 
Factor 





Figure IV-4 ELF 86 optical field, electron beam diagram 


The ratio of the electron beam area to the optical beam 
area decreased slightly giving a filling factor F=0.024, as 
seen in Figure IV-4. The electron pulse length remained at 


20 slippage distances with minimal short pulse effects. The 
dimensionless beam size o, remained nearly the same due to a 


slightly smaller optical wavelength 0.72 cm: The 
dimensionless current increased to j=5000. The input optical 


field is provided by a 100kW, 35GHz magnetron giving a,=14. 
The emittance was significantly reduced to £-0.07 cm- 


rad. The overall increase in current and better beam quality 


has increase the power output from 80 MW for ELF85 to 180MW 


for ELF86 and the efficiency increased from n=5% to n=6%. 


G. ELF87 
A major change in ELF87 was the use of a high-current 


induction linac. This linac provided a beam energy of 3.5 
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MeV and a peak current of 1 kA. The dimensionless current 


1000 and the undulator was also extended to L=4 m. 


HE 


Another major change was the shape of the electron 


pulse. The beam length cross section was oval as seen 
Figure IV-5. However, this beam shape has little effect 


the model. 
The shift in the optical wavelength to 7ай 175 


causes the optical mode to become much smaller. The ratio 


electron beam area to optical mode area is much smaller. 


in 


on 


cm 


ФЕ 


А 


process called optical guiding also causes the optical mode 


to remain small. 





Figure IV-5ELF 87 optical field, electron beam diagram 


The net effect is a larger filling factor, F=0.24, and more 


interaction between the light and the electron pulse. 


Не MIT FEL 


The MIT FEL is operated by the Department of Physics 


and Research Laboratory of Electronics, Massachusetts 


SL 


Institute of Technology. ThiserEL was deésitoned @temererare at 
millimeter and sub-millimeter wavelengths. The MIT FEL has 
not been studied as closely, since it has lower gain than 
ELF. The most significant difference between this laser and 
the ELF laser is the undulator. The MIT laser uses a helical 
undulator. 


L- Natural Focusing 


An advantage of the helical undulator is the natural 
focusing of the electron beam. The magnetic field is always 
pushing the electrons toward the center of the beam. Far an 
x-ray FEL, any added benefit from natural focusing could 
result in a significant increase in gain. 

2 SASE 

The other significant feature of the MIT laser is the 
SASE startup at long wavelengths. From its operation, it can 
be seen that SASE occurs within about the first 20 undulator 
periods. After the first 20 periods, the FEL field growth 
becomes exponential and can be simulated using a small 
initialweptical field да. | 

cm Electron Beam 

A multielectron gun with a Pulserad 110A accelerator 
creates the electron beam. The beam is of relatively low 
kinetic energy, 2.3 MeV, with a peak current of 930 amps. 


The pulse length is 7200 times the slippage distance, so 
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БЕС smEEpulsecmescsedebs wath tha electrons along the 


entire length of the undulator. 


Ts CONCLUSION 

In careful study of these FEL's, one can discern which 
parameters have the most influence on the growth rate and 
FEL interaction. A simulation model can often be simplified 
by concentrating on those parameters which create the large 
growth rate or decays, and ignoring those which have no 
discernible effect on the electron/optical field 
interaction. The parameters developed in this chapter can be 


used in a simple model. This is left for further study. 


Ба 
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V. DARMSTADT FEL 


The next section is an analysis of the Darmstadt FEL 


utilizing cCömputer “Simulations. The Darmstadt FEL is a 
medium gain laser and operates between A=6.6 and 7.8um 


wavelength, with electron beam energies between 32 and 29 
MeV. It is driven by the S-DALINAC, а re- circulating. 
electron linac using superconducting acceleration cavities, 


which allow for a continuous wave electron beam. 


A: DARMSTADT PARAMETERS 


The parameters for the Darmstadt FEL are shown in Table 
V-1. The primary use for the Darmstadt FEL is to provide 
infrared light for medical research. Most of- the 
experiments involve medical uses such as ablation of 
biological tissue. 


Electron Energy-KE 
Peak Current-I 
Electron Pulse Length-l, 
Number of Periods-N 
Undulator period-A, 
Undulator Parameter-K 
Dimensionless current-j 


Energy Spread-AE/E 
Q-factor 





Table V-1 Darmstadt parameters 
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B. BEAM SIZE EFFECTS ON GAIN 


Initially the parameter values obtained for the 
Darmstadt FEL were not properly understood. The electron 
beam radius was given at the end of the undulator, and it 
was mistaken to be the radius at the center of the 
undulator. After several communications and some analysis, 
the discrepancy was discovered. As a result of angular 
spreading along the undulator the beam radius was 0.2mm at 
the waist, and a little less than imm at the end of the 
undulator. For an electron beam with Darmstadt's parameters 


and an electron beam waist radius r,-1mm, the angular 
du er sole ion would Бе 657-227 Іп “ре time it "took ЕО 


discover the error, an analysis of the influence of angular 
spread on gain was conducted. Despite not having апу 
relevance to the Darmstadt FEL, the analysis was interesting 
and may have future relevance. Therefore the results are 
included here. The analysis below was conducted assuming an 
electron beam waist radius of 1mm. 


Je Theoretical Background 


An angular distribution әш бе) 5 is signi icant enough 


to degrade the gain of an FEL. A large electron beam radius 
means that there are a significant number of the electrons 
which are injected into the undulator that are off-set from 


the axis of propagation and at an angle with respect to that 
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axis. The larger the diameter, the more electrons that are 
injected off center and at larger angles. The electrons on 
the outer portion of the beam experience a larger magnetic 
field. These outer electrons are "wiggled" harder апа 
therefore travel on a slightly longer path. As a result of 
the longer path, the electron slows down in the longitudinal 
direction with respect to an electron in the center of the 
beam. The spread in longitudinal velocities creates a spread 
in phase velocity. A spread in phase velocity then results 
in a redugegeon in gain: | 

The same argument can be used for electrons that are 
injected at an angle with the axis. An electron injected at 
an angle moves toward (or away from) the more intense 
magnetic field off axis. The electron is pushed back to the 
center harder than the other electrons and travels a longer 
distance. The longer distance means a greater spread in 
phase velocity and less gain. 


2 Method of Analysis 


In order to determine the effect of beam size on the 
gain, the effect of value of og on the gain is analyzed. The 
gain for each og must be calculated at the value of 


desynchronism giving maximum gain. For each og, there is a 


different desynchronism value d for which gain is a maximum. 


The amount of time to run enough simulations to determine 


the maximum gains for each og would have been prohibitive. 


57 


Gain versus desynchronism curves were plotted for 09=7.3 and 
Go-2, Nas (Shown in Figure SUE ІП “order Jee, tee this 


information to determine how the maximum d depends on Gg, 
linear interpolation was used. Assume a linear relationship 


between the maximum desynchronism d and oe using the form 
Jus c tora sreapsgie mmm gain for each 
value of og can be estimated in subsequent simulations. 


Maximum gain was found for 09=7.3 at d=0.003, and for 0$g23 


at d=0.0025. After a small amount of algebra we find that 


m > 0200058. таме" 520200682 Io сезе спе алау сйс 
estimation, Og=5 was estimated to be a maximum at d=0.0037. 


A series of simulations shows that in fact, the maximum 
occurred at 0.0036; in good agreement with the interpolated 
formula. In Figure V-1 it can be seen that there is a flat 
plateau about the maximum of the curve, so our c SEE 


for the maximum is accurate enough. In 

Figure V-2 the gain found for each oe, at the maximum d 
determined by the interpolation formula above, is plotted. A 
value of 0g=0 means perfect injection, no deviation іп 
either angle or position. The actual beam radius r,=0.2mm 


corresponds to o9=0.3 which is nearly perfect. The gain for 


0920.5 is about 3% less than for perfect injection which is 


an insignificant reduction. 
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Figure V-1 Gain Versus Desynchronism 


On the other hand, the original value of 0e-7.3 results in a 


77% decrease in gain which is not negligible. In fact, it 
would mean the Darmstadt FEL would not work. 

From Figure V-2 the effects of an angular or radial 
spread on the reduction of Darmstadt gain can be straight- 


forwardly determined. 
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e 


Gain 





Og 


Figure V-2 og versus Gain 


The final power P versus desynchronism d curves for 09=7.3 


and 09-3 were plotted, as seen in Figure V-3. The final 


power in strong fields at saturation peaks at small values 
of d where gain is small, and extends out to larger values 


E ci 
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Power 


0 0.002 0.004 0.006 0.008 0.01 0012 0.014 0.016 02078 
d 


Figure V-3 Power Versus Desynchronism 


er TAPER EFFECTS 


Using the actual Darmstadt FEL electron beam radius of 
7. -0220m, өле “спе "Грататесете топ Table vV- I тептер г 
versus desynchronism graphs for several tapers were 


determined. The following positive and negative tapers were 


analyzed: 6--5.4л; -2.5Л, Össi ns 47 


Ls Compiled Data 


Figure V-4 and Figure V-5 are 3-D graphs that depict 


power (on the vertical axis) as a function of both taper 6 


ба! 





0 


20/2 


enl zh + > 3 = —À S f? die e y сж 4 ^ Sa RT 
ana desyncnhnron-sm с. Tre two figures contain the same data 


but presented irom gdirceren- points of view. 





Figure V-4 3-D Plot of Power Versus, d, and 6(Positive View) 

The power versus c curves P(a) for each value of taper 
ð have similar structures. When 8-0, there is no power 
output at any 6. At approximately a=0.008, che output power 
is maximum for each taper 6. The cesynchronisms in tne range 


0<4<0.008 dc not prosuce any power for any taper. 


с) 





Figure V-5 3-D Plot of Power, d, and ó(Negative View) 


The sharp rise in power and high peak power is typical of 


many FELS cperating in this parameter regime. Tne largest 


peak power occurs at 6=0. 


A high quality factor Q means There are small optical 


Loses per pass throuch the undulator. These losses include 


mirror loss and mirror transparency wnere optical 
removed for tne end user. A high Q leads to sharper 
in the cesynchronism curves ana to other important 


ч 


explained ir the next section. 
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power is 


22 Untapered Power Curve Analysis 


To show more detail of the data in Figure V-4 and 


Figure V-5, Figure V-6 and Figure V-9 show power versus 


desynchronism with no taper 6=0 and the highest value of 


negative taper д--5.4л. 


600. 


200 


300 


Power 





200 


100 





0 0005 00 0.015 0.02 0-025 Ооз 02035 0.04 


Figure V-6 P Versus Desynchronism 0-0 
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ПОЕТЕСИ СЕС ЕЕ exhibrtsesuteevele behavior 
for desynchronisms between d=0.0008 and d=0.011 (ref (2)), 
and thus is not stable. Since the power levels are not 
stable the power in this region is an average of the power 
oscillations. 

m Simulation 

Figure V-7 is the output from a single simulation run 
for the Darmstadt FEL. The input variables for this run are 
displayed in the top square. There are nine output windows 
displayed below the input window. 


The top-left window shows the dimensionless optical 
field amplitude latz) | at the end of n=6000 passes. The 


middle-left window shows the evolution of the field 





amplitude after at each pass la(z,n) , in this case for 


n=6000 passes. The lighter the color, the higher the 


amplitude, as is the case for each middle window. In the 
bottom-left window, the parabolic electron pulse j(z-t) is 


plotted for reference at the beginning and end of the 
undulator; their separation shows the slippage along the 
meum aror. 


The bottom-center window plots for reference the weak- 
field gain spectrum as a function of phase velocity, v. The 
middle-center window shows the optical power evolution 


P(v,n) as a function of phase velocity v at each pass n. The 
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jindow show the final opticai power spectrum Pí(v) 


s 


top-center 
at tne end of n=6000 passes. 

The pnocttom-right wincow shows the total optical power 
evolution D(a) as a function of the pass mumber n. The 


center-right wincow shows the evolution of the electron 


phase velocity and number of passes n. The upper-right 
window displays che final electrcnü phase ewe 


Poem OUL Ome (VP oO00) ae ele ecqor neo000 Dasses. 


О, 


| **** FEL Pulse Evolution **** 
'3=1.6 6,-1.1 6,40 6,=2.3 d=0.006 
loz1t230 A= А=0 б=0 N-80 8{=0.0001 


P(v,n) | е, Пп) 


teen ma аласа 2 сми 





Figure V-7 Trapped Particle Instability 


4. Limit-Cycle Behavior in the Untapered case 


In the  untapered  undulator, the trapped-particle 
instability occurs for small values of desynchronism at 
d=0.012. In strong optical fields, electrons become trapped 
in the deep potential wells in phase space. These trapped 
electrons oscillate and radiate into side- bands off the 
main wavelength of the optical field. In the center window 
of Figure V-7, the side-bands off the main peak can be seen. 
Due to the power in the side-bands, limit-cycle behavior is 
present in the total power P(n), and can be seen in the 
output power shown in the lower right hand window. The 
output power oscillates between two power levels as the FEL 
oscillates between two frequencies. The trapped particle 
instability does not always cause limit cycle behavior. In 


Figure V-8, the power exhibits chaotic output. 
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| **** FEL Pulse Evolution **** 
1-1.6 6.51.1 69-0 Og=2.3 d-0.004 
| Q-110 A=0 A=0 0-0 N-80 06.20.0001 


P (v,n) 
r1! 












0.21? (п) 297.1. 
| + 


t 
it 





n 10000 


Figure V-8 Chaotic Power Output 


SE Limit-Cycle Behavior in the Negative Tapered Case 
cure V-9 is the desynchronism curve P(d) for the case 


or this taper O-2-5.4nz, limit 


т] 


of negative taper, ð=-5.47. 
cycle behavior never occurred, or was so smali that it was 


not noticeable. 
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Figure V-9 P versus d 6=-5.4 


In Figure V-10, at the specific desynchronism d=0.005 
and negative taper 6=-5.4 we see there are no trapped 


particles. Above a certain value for d, d=0.004, the trapped 
particle instabilities do not occur. This results in smooth 
power output with no side bands around the optical 
wavelength as seen in the center window. The overall effect 
is a lower power output overall but provides for cleaner 
output power. This may be important for some applications 
where a narrow spectrum is important. For other applications 


more power may be important. 
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**** FET Pulse Evolution **** 
'j=1.6 0,=1.1 09=0 Og=2.3 d=0.005 
Q=110 =0 =0 9=-5.47 №80 96=0.0001 
Ғ (У, п) 
i Y 


Figure V-i0 No Trapped Particles 


In Figure V-11, the efiects of the trappec particle 
instability can be seen in the left and center windows. 
There is no limit cycle behavior as is seen from the steady 
power P(n) window in the bottom right-hand window. In every 
caper studied, there were values of а for which the trapped 
particle :mctabsxity occurred Опе no Linit cycies. 


жжжж 





FEL Pulse Evolution **** 
(351.6 6,-1.1 660 0с=2.3 4-0.003 
Q-110 А-0 A=0 


o=-5 . 44 N=80 


fi 55.48 





T" 
Figure V-11 Trapped Particles with Stable Power 


D: CONCLUSION 

The data collecteä £rom simulations of the Darmstadt 
FEL describe how Gecreasing electron Deam quality would 
decrease that 


observed gain. Also simulations showed that 
there are 


Operation 


ana chaos may be expectec 


regions of where limit-cycle behavior 
tO occur in the experiment. If was 
also observed that the i 


i 


( 


rapped-particie instability may 
Occer without resence of limi 


it-cycles. 
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